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The observation of new suppressed B° s decays, B® — > J/^fK* (892)° and B° s — > J/wK®, and the 
measurement of their branching ratios is presented. This measurement is based on an integrated 
luminosity of 5.9 fb _1 of CDF data collected by a dedicated di-muon trigger. A cut based opti- 
mization is carried out for the observation of B® — > J/\j/K*°, while a neural network is used for 
the B® — s- J/y/K®. In addition to the observation of the new decay modes, the ratios of branching 
fractions to the reference B° decays are measured. 
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While B decays have been extensively studied at the B factories experiments, much less is 
known about B s decays. This contribution presents studies of two specific B s decays [jl|]: B s — > 
J/yK® and B° s J/yK* (892)° , with J/y -> pi + H~, K% -> k+k' and K*(S92)° -> K %. -> 
J/yK® is a CP eigenstate and has never been observed. Measurement of its lifetime directly probes 
the lifetime of the heavy mass eigenstate, %b s h ■ Additionally, large samples of B® — > J/yK® can 
be used to extract the angle 7 of the unitary triangle The B° — > J/yK decay is yet another 
unobserved mode which contains an admixture of CP final states. With a larger data sample, an 
angular analysis of B® — > J/yK*® can be carried out to help disentangle penguin contributions in 
B® — > J/y<p [§]. In addition to the first observation of these two decays, the ratios of branching 
ratios of B® — > J/yK® and B®, — > J/yK*® to the reference B° decays are measured using the relation 

SS(B® -> J/yK)/m(B° J/yK) =A rel x f d /f s x N(B° S -> J/yK)/N(B° -> //yrtt), 

where ^ represents ^Tj or By measuring the ratio of the number of decays, N(B® — > J/yK) 
and N(B° — >• J/yK), from data and the relative acceptance, A re /, between the fi° and B® from Monte 
Carlo simulation (MC), the value 3${B° S -> A") /^(fi° -)■ is extracted by inputting the 

ratio of fragmentation fractions f s /fd- 

The data used in these analyses are selected from a sample enriched in J/y — > }JL + pL~ decays, 
collected by the CDF Run II detector [ffl]. The integrated luminosity of this sample is 5.9 fb _1 . 
The J/y dataset contains events with at least one reconstructed J/y selected by dedicated di-muon 
triggers. In addition to the selected J/y, two tracks are combined with them in a kinematic fit to 
reconstruct 5° ->■ J/yK® and B° ->■ J/y K*° candidates. For the B° ->■ J/yK% analysis, the two 
tracks are reconstructed as pions and combined to define a K$ candidate. The K* candidate for the 
B° — > J/yK*° decay is reconstructed from the combination of a % and a K. 

The final event selection in the B° — > J/yK*° analysis is optimized by maximizing 5/(1.5+ VB). 
This quantity is well suited for signal discovery as described in A simultaneous four-dimensional 
optimization is carried out over 4 quantities: Pt(k), Pt(K), transverse decay length 7^ (6°) andB^ 
vertex fit probability. For the purpose of extracting the yields of B° — > J/yK*° and B° — > J/yK* Q 
signals in the invariant mass distribution, an accurate modeling of signals and backgrounds is 
needed prior to the fit. The signal contributions are modeled with three Gaussians template ob- 
tained from a fit to B° MC. The B*' template used in the final fit is identical to B° template, except 
for a shift of 86.8 MeV/c 2 in the mean value of the three Gaussians. This value corresponds to the 
known value ^ for the mass difference between B® and B°, A(m B o — m B o). The backgrounds con- 
sidered in this analysis are combinatorial background, partially reconstructed b hadron decays and 
6° — > J/y<p decay. The first one is resulting from different sources, for example a real J/y plus two 
random tracks, where the J/y could be a prompt J/y or coming from a B decay. Other sources that 
could contribute to it are fake J/y reconstructed with prompt fake muons or fake muons coming 
from heavy flavor. The combinatorial background is modeled in the final fit with an exponential 
function. The partially reconstructed background, which is fitted with an ARGUS function [Rj, is 
partially reconstructed b hadrons where a five-body decay occurs where a %, K, or 7 is not recon- 
structed. Finally, to model the B° — > J/y<p background, a template consisting of two Gaussians, 
extracted from simulation, is used. The size of B® — > J/y<p contribution is constrained using data. A 
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binned log likelihood fit is performed to the invariant mass distributions using the templates for sig- 
nals and the functions described above. Figure [j] shows the fit, including the different contributions. 
From the fit, the yields of the B° -> J/yK*° and 5° -> J/yK*° signal are 9530 ± 1 10 and 151 ± 25, 
respectively. The statistical significance of the B® — > J/yK*° signal is 8.0a. The measured ratio 
of yields, N(B° S -»• J/yK*°)/N(B Q -»• J/yK*°), is 0.0159 ± 0.0022 (stat.)±0.0050 (sys.). The sys- 
tematic uncertainty is dominated by the combinatorial background contribution uncertainty, with 
a relative uncertainty for the ratio of 31.4%. The other sources of systematic uncertainty are the 
signal modeling (4.4%), uncertainty on A(m B o —m B o) (0.1%), combinatorial background modeling 
(1.3%) and B° s ->■ J/Y<t> contribution (1.3%). 
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Figure 1: Invariant mass distribution for J/y/K* and fit including the different contributions (left). The 
disttibution is enlarged in the signal region for more detail (right). 

The B° — > J/\yK® decay has two main differences with respect to the B° — > J/yK* decay. 
First, it contains a which has a long lifetime. The displacement between the K% candidate 
vertex and the B candidate vertex is used in the event selection to reduce backgrounds from prompt 
sources. Moreover, the K® is narrow resonance, so a mass constraint can be applied and improve 
mass resolution. Second, the B®, — > J/yK® signal contribution is expected to be smaller than the 
B® — > J/yK*° signal contribution. Therefore it is crucial to effectively suppressing combinatorial. 
A neural network (NN) is used a discriminator. In order to train the NN, simulated B® events are 
used as signal. Data from the upper side band in the B° invariant mass distribution, well separated 
from the signal region, are used as a background data sample. Twenty-two different kinematic 
and topological variables are chosen as inputs for the NN training. The selection is optimized 
by maximizing 5/(1.5 + y/B) and a likelihood fit is performed to the invariant mass distribution 
to extract the yields of B° — > J/yK® and B° — > J/yK® signals. Signals and the two backgrounds 
are modeled with the same functional form as in the B° — > J/\yK*° analysis. Moreover, in the 
B s — > J/yK® analysis , At, — > J/\j/A, where A — > p%, is a background when the p is reconstructed 
as a n. This contribution is suppressed using a cut in the angular distribution between K% and the 
K with lower p T . From the fit, shown in Figure |2[ the yields of the B° — > J/yK® and B® — > J/yK^ 
signal are determined to be 5954 ± 79 and 64 ± 14, respectively. The statistical significance of the 
5° J/yK° s signal is 7.2a. The value of N(B® J/yK°)/N(B° -> J/yK° s ) is 0.0108 ± 0.0019 
(stat.) ± 0.0010 (sys.). The sources of the systematic uncertainty are similar to the other analysis. 
In this case the relative uncertainties for the ratio are 5.6% from the combinatorial background 
contribution, 5.6% from the combinatorial background modeling, 4.6% from the signal modeling 
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and 0.1% from the A(ra B o — ra g o). 
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Figure 2: Invariant mass distribution for J/yK^ and fit including the different contributions (left). The 
distribution is enlarged in the signal region for more detail (right). 

To determine the 3S(B° S -> J/yK) /3$(B Q J/yK), where K represents K§ or K*°, the relative 
acceptances of B° -)■ J/yK$ to 5° -> 7/^° and B° -> //V^T to 5° -> J/yK*° need to be deter- 
mined from MC. The value for A re i is determined to be A ret = 1.012 ± 0.010 (stat.) ± 0.042 (sys.) 
for the K$ analysis and A rd = 1.057 ± 0.010 (stat.) ± 0.263 (sys.) for the A"* analysis. The £° 
and fi° lifetimes, B hadron pj spectrum and polarization, this last one only for the B° s — > J/\j/K* Q 
analysis, induce uncertainties in A re i. 

To determine f s / fd > the most recent CDF measurement [§] of f s / (/„ +fd) x 88(D S —?■ <pn) is 
combined with the actual known value for 3$(D S — > (fin). With the input of f s /fd = 0.269 ± 
0.033, the ratio of branching fractions to the reference B° decays are: 

m{B Q s -> J/yK*°)/^(B° -> = 0.062 ± 0.009(tfaf.) ± 0.025(jyj.) ±0.008(/rag.), 

-> J/yKl)/£g(B Q = 0.041 ±0.007 (jfitf.)± 0.004(.syj.) ± 0.005 (/rag.). 

The world-average values for ^(B — >■ J/yK*°) and ^(2?° — >• J/yK°) are used for normalization 

to calculate the absolute branching fractions: 

m(B° s -> J/\ifK*°) = (S.3±\.2(stat.)±3.3(sys.)±l.0(frag.)±0A(Norm.)) x 10~ 5 , 
£g(B° ^-J/yK ) = (3.5 ±0.6(stat.)±0A(sys.)±0A(f rag.) ±0A(Norm.)) x 10~ 5 . 

This measurement is yet another CDF contribution to the exploration of bottom-strange mesons. 
Further are expected with the 10 fb -1 data sample expected in a year from now. 
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